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Abstract 

Following the prescription of Ref. [1] in which perfect state transference (PST) of a 
qubit over distance regular spin networks was discussed, in this paper PST of an arbi- 
trary d-level quantum state (qudit) over antipodes of more general networks called pseudo 
distance-regular networks, is investigated. In fact, the spectral analysis techniques used 
in the previous work [1], and algebraic structures of pseudo distance-regular graphs are 
employed to give an explicit formula for suitable coupling constants in the Hamiltonians 
so that the state of a particular qudit initially encoded on one site will evolve freely to 
the opposite site without any dynamical control, i.e., we show that how to derive the 
parameters of the system so that PST can be achieved. 
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The transference of quantum information, encoded in a quantum state, from one part of a phys- 
ical unit, e.g., a qubit, to another part is a crucial ingredient for many quantum information 
processing protocols [2]. There are various physical systems that can serve as quantum chan- 
nels, one of them being a quantum spin system. Quantum communication over short distances 
through a spin chain, in which adjacent qubits are coupled by equal strength has been studied 
in detail, and an expression for the fidelity of quantum state transfer has been obtained [3, 4]. 
Similarly, in Ref. [5], near perfect state transfer was achieved for uniform couplings provided 
a spatially varying magnetic field was introduced. After the work of Bose [3], in which the 
potentialities of the so-called spin chains have been shown, several strategies were proposed to 
increase the transmission fidelity [6] and even to achieve, under appropriate conditions, perfect 
state transfer [7, 8, 9, 10, 11, 12]. All of these proposals refer to ideal spin chains in which 
only nearest-neighbor couplings are present. In Refs. [7, 8], the d-dimensional hypercube with 
2*^ vertices has been projected to a linear chain with d + 1 sites so that, by considering fixed 
but different couphngs between the qubits assigned to the sites, the PST can be achieved over 
arbitrarily long distances in the chain. In Ref. [1], the so called distance-regular graphs have 
been considered as spin networks (in the sense that with each vertex of a distance-regular 
graph a qubit or a spin was associated ) and perfect state transference (PST) over them has 
been investigated, where a procedure for finding suitable coupling constants in some particular 
spin Hamiltonians has been given so that perfect transference of a quantum state between 
antipodes of the networks can be achieved. The aim of this paper is to extend this proposal to 
systems of particles with arbitrary number of levels, the so-called qudits. These systems can be 
appeared in condensed matter and solid state physics such as the fermionic SU (N) Hubbard 
model [13, 14, 15]. In Ref. [16], state transference over spin chains of arbitrary spin has been 
discussed so that an arbitrary unknown qudit be transferred through a chain with rather good 
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fidelity by the natural dynamics of the chain. In this work, we focus on the situation in which 
state transference is perfect, i.e., the fidelity is unity. Furthermore, we consider more general 
graphs called pseudo distance regular graphs or QD type graphs [17, 18] (distance regular 
graphs are special kinds of pseudo distance regular ones) as underlying networks and investi- 
gate PST over antipodes of these networks. In fact this work is an extension of the previous 
work [1] to PST of a qudit over a pseudo distance regular network. We use the techniques 
such as stratification of the graphs [17], [19]- [23] and spectral distribution associated with the 
graphs. Then, we consider particular hamiltonians with nonlinear terms and give a method 
for finding a suitable set of coupling constants so that PST over antipodes of the networks 
be possible. In fact, we give an explicit formula for suitable coupling constants so that PST 
between the first node of the networks and the opposite one can be achieved. As examples we 
will consider some important pseudo distance regular networks such as the G„ networks, the 
modified G„ networks, Hadamard network, etc. 

The organization of the paper is as follows: In section 2, we review some preliminary facts 
about graphs and their stratification, pseudo distance-regular graphs and spectral distribution 
associated with graphs. Section 3 is devoted to PST of a qudit over antipodes of pseudo 
distance-regular networks, where a method for finding suitable coupling constants in particular 
spin Hamiltonians so that PST be possible, is given. The paper is ended with a brief conclusion 
and one appendix. 

2 Preliminaries 

In this section we recall some preliminaries related to graphs, their stratifications and the 
notion of pseudo-distance-regularity (as a generalization of distance regularity) of graphs. 
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2.1 Graphs and their stratifications 

In this section, we review the stratification of the graphs and the notion of pseudo-distance- 
regularity. 

A graph is a pair F = (V, E), where F is a non-empty set called the vertex set and E is 
a subset of {{x,y) : x,y E V,x y} called the edge set of the graph. Two vertices x,y E V 
are called adjacent if (x, y) e E, and in that case we write x ~ y. For a graph F = (V, E), the 
adjacency matrix A is defined as 

r 1 if a - /3 

iA)a,fS = . (2-1) 

otherwise 

Conversely, for a non-empty set V, a graph structure is uniquely determined by such a matrix 
indexed by V. 

The degree or valency of a vertex x eV is defined by 

K{x)^\{yeV:yr^x}\ (2-2) 

where, | • | denotes the cardinahty. The graph is called regular if the degree of all of the vertices 
be the same. In this paper, we will assume that graphs under discussion are regular. A finite 
sequence Xo,Xi, ...,Xn G ^ is called a walk of length n (or of n steps) if ~ Xi for all 
i = 1, 2, Let f{V) denote the Hilbert space of C-valued square-summable functions on 

V. With each (3 E V we associate a vector such that the /3-th entry of it is 1 and all of 
the other entries of it are zero. Then {\P) : P E V} becomes a complete orthonormal basis of 
l'^{V). The adjacency matrix is considered as an operator acting in P{V) in such a way that 

A\/3) = Y: (2-3) 

Now, we recall the notion of stratification for a given graph F. To this end, let d{x,y) 
be the length of the shortest walk connecting x and y ior x ^ y. By definition d{x, x) = 
for all X E V. The graph becomes a metric space with the distance function d. Note that 
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d{x, = 1 if and only if x ~ We fix a vertex o e ^ as an origin of the graph, called the 
reference vertex. Then, the graph F is stratified into a disjoint union of strata (with respect 
to the reference vertex o) as 

oo 

V=[j r,(o), r,(o) -.^{aeV: d{a, o) = i} (2-4) 

i=0 

Note that Ti(o) — may occur for some i > 1. In that case we have Ti(o) — rj+i(o) = ... = 0. 
With each stratum Ti{o) we associate a unit vector in 1^{V) defined by 

1'?'^) = ^ E I")' (2-5) 
where, = |rj(o)| is called the i-th valency of the graph (/tj := |{7 : (9(o, 7) = i}\ = \Ti{o)\). 



2.2 Pseudo-distance regular graphs 

Given a vertex a e of a graph F, consider the stratification (2-4) with respect to a such that 
ri{a) = for i > d. Then, we say that F is pseudo-distance- regular around vertex a whenever 
for any /3 G Tk{a) and < k < d, the numbers 

(2-6) 

do not depend on the considered vertex 13 e Ffe(Q;), but only on the value of k. In such a case, 
we denote them by c^, and hk respectively. 

It should be noticed that for regular graphs F (k(/3) — k = Ki for all (3 &V), the numbers 
Cfe, afe and hk read as 

Cfc=|Fi(/3)nFfe_i(Q;)|, afc = |Fi(/3)nFfc(Q;)|, 6^ = |Fi(/3) n Ffe+i(Q;)|, (2-7) 

where we tacitly understand that F_i(q;) = 0. The intersection numbers (2-7) and the valencies 
Ki = |rj(a)| satisfy the following obvious conditions 

(^i bi -\- Ci = K, Ki—ibi—i = ACjCj, i = 1, (i. 
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Ko = ci = 1, bQ = Ki = K, (co = bd = 0). (2-8) 

One should notice that, the definition of pseudo-distance regular graphs together with the 
Eq.(2-8), imply that in general, the valencies Ki ( the size of the i-th stratum) for i = 0, 1, d 
do not depend on the considered vertex (3 e I'kicx), but only on the value of k. 

The notion of pseudo-distance regularity has a close relation with the concept of QD type 
graphs introduced by Obata [17], such that for the adjacency matrices of this type of graphs, 
one can obtain a quantum decomposition associated with the stratification (2-4) as 

A^A+ + A- + A^, (2-9) 

where, the matrices A'^, A~ and ^4° are defined as follows: for /? e rk{a), we set 

^+1/?)= y: e ^°i/5)= e (2-10) 

<5Grfe+i(a),5~/3 5erfe_i(a),<5~/3 SeTk{a),S^l3 

Since /3 e ^k{(^) and /3 ~ 5 then S e ri._i{a) U Tk{a) U rk^i{a). 

It has been shown in Ref. [17] that, if the closed subspace of l'^{V) spanned by = 
0, . . . , d — 1} be invariant under the quantum components A'^, A" and then there exist two 
sequences (called Szego- Jacobi sequences) {a;;}^^ and {q;;}^^ derived from A such that 

A-|0o) = O, A-\<Pi) = ^i\<Pi_i), />1, 

A''\(t>i)^ai\<t>i),l>Q, (2-11) 

where ui+i = ^k?_(j) with K_(j) = |{i e r/(Q;) : i ~ j}l i ^ r;+i(Q;), and a; = «;o(j) = 
|{i e VJ;i ~ j}| for j e r;(Q;), for I > 0. One can easily check that the coefficients ai and loi 
are given by 

ai^ K-bi- ci, ui = Z = 1, d. (2-12) 

By using (2-9) and (2-11), one can obtain 

= + + ^>o. (2-13) 
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One should notice that, the vectors \4>i),i = 0,1, ...,d— 1 form an orthonormal basis for the so 
called Krylov subspace Ka{\(l)o),A) defined as 

Xd(|0o),A) = span{|0o),A|0o),---,^'~'|0o)}. (2-14) 
Then it can be shown that [24], the orthonormal basis are written as 

|0,) = P,(A)|0o), (2-15) 

where Pi — aQ + aiA + ... + aiA'^ is a polynomial of degree i in indeterminate A (for more details 
see for example [19, 24]). 

It may be noted that, the pseudo-distance- regularity is a generalization of the notion of 
distance-regularity which is defined as: 

Definition (distance-regular graphs). A pseudo-distance regular graph F = {V, E) is 
called distance-regular with diameter d if for all k e {0, 1, d}, and a, (5 & V with e ^k{(^), 
the numbers Ck{l3), ak{P) and bk{P) defined in (2-6) depend only on k but do not depend on 

the choice of a and p. 

Now, it should be noticed that, the stratification of distance- regular graphs will be indepen- 
dent of the choice of the reference vertex (the vertex which stratification is done with respect 
to it). 

2.3 Spectral distribution of the graphs 

It is well known that, for any pair (^4, \(f)o)) of a matrix A and a vector \(f)o), one can assign a 
measure /i as follows 

Mx) = (</.o|£^(x)|0o), (2-16) 

where E{x) = J2i \ui){ui\ is the operator of projection onto the eigenspace of A corresponding 
to eigenvalue x, i.e., 

A= [ xE{x)dx. (2-17) 
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Then, for any polynomial P{A) we have 

P{A) = J P{x)E{x)dx, (2-18) 

where for discrete spectrum the above integrals are replaced by summation. Therefore, using 

the relations (2-16) and (2-18), the expectation value of powers of adjacency matrix A over 
reference vector |0o) can be written as 

{(po\A"'\<Po) = ! x^ti{dx), m = 0, 1, 2, .... (2-19) 

JR 

Obviously, the relation (2-19) implies an isomorphism from the Hilbert space of the stratifi- 
cation onto the closed linear span of the orthogonal polynomials with respect to the measure 
//. 

Prom orthonormality of the unit vectors given in Eq.(2-5) (with |0o) as unit vector 
assigned to the reference node) we have 

% = {(t>i\h) = I Pi{x)P,{x)fi{dx), (2-20) 

where, we have used the equation (2-15). Now, by substituting (2-15) in (2-13) and rescaling 
Pk as Qk = s/uJiTTTuJkPk-i the spectral distribution fi under question will be characterized by 
the property of orthonormal polynomials {Qk\ defined recurrently by 

Qo{x) = 1, Qi{x) = X, 

xQk{x) = Qk+i{x) + akQk{x) + uJkQk-i{x), k>l. (2-21) 

If such a spectral distribution is unique, the spectral distribution /i is determined by the 
identity 

iRX-y x-ao- ^2 Qd+i{x) f^n x - xi 



x—a\- 

x—oto- 



where, xi are the roots of the polynomial Qd+\{x)- Gp}yx) is called the Stieltjes/Hilbert trans- 
form of spectral distribution [x and polynomials {QI^''} are defined recurrently as 

Qf{x) = \, Qf{x)^x-a^, 
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xQ^\x) = ggi(x) + ak+,Qi\x) + Uk+iQ^Ux), k>l, (2-23) 
respectively. The coefficients ji appearing in (2-22) are calculated as 

7/ [(^ - xi)G,{x)] (2-24) 

Now let G^{z) is known, then the spectral distribution fi can be determined in terms of 
X/, Z = 1, 2, ... and Gauss quadrature constants 7;, 1 — 1,2, ... as 

d 

l^ = J2^iS{x-xi) (2-25) 

1=0 

(for more details see Refs. [18, 25, 26, 27]). 

3 PST of a qudit over antipodes of pseudo distance- 
regular networks 

3.1 State Transference in c?-dimensional Quantum Systems 

A d-dimensional quantum system associated with a simple, connected, finite graph G — (V, E) 
is defined by attaching a ci-level particle to each vertex of the graph so that with each vertex 
i eV one can associate a Hilbert space Hi ^ C^. The Hilbert space associated with G is then 
given by 

Hg = ^^.yUr = {CT'', (3-26) 

where := \V\ denotes the total number of vertices (sites) in G. 

Then, the quantum state transfer protocol involves two steps: initialization and evolution. 
First, a quantum state 

d-l 

\i>)A = ao\0)A + a^W)A e Ti-A 

v=l 
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(with e C and Ylt=o l^i^P = 1) to be transmitted is created. The state of the entire spin 
system after this step is given by 

\jj;(t = 0)) = |V'a)(8)|0...00b) = ao\OA)®\O..mB)+ai\lAMO..mB)+...+ad-i\(d-l)A)0\O..mB). 

(3-27) 

Then, the network couphngs are switched on and the whole system is allowed to evolve under 
U{t) — e~*^* for a fixed time interval, say to- Prom the fact that H\Oa) <S> |0...00b) = 0, the 
final state at time to will be 

d-l N 

mto)) = ao|OAO...OOB) + E «4E /fc^(^o)|0 . . . ^ 0...0)}, (3-28) 

^^=1 ^==1 k-th 

where ftlito) (0...0^0...0|e-'-^*o|z/AO...O) for k = 1,2, ...,7V; = 1, . . . , d - 1. In order 

k-th 

to transfer the state IV'a) to the site B perfectly (in order to PST is achieved), the following 
conditions must be fulfilled 

l/is(^o)| = 1 for i/ = 1, 2, d-l and some < to < oo (3-29) 

which can be interpreted as the signature of perfect communication (or PST) between A and 
B in time to. The effect of the modulus in (3-29) is that the state (3-28) will be 

d-l 

IV'(to)) = ao|OAO...OB) + Y,e''^''au\OAO...O) ® \iy)B 
so, the state at S, after transmission, will no longer be |V')a, but will be of the form 

The phase factors e*"^" for v — 1, 2, d — 1 are independent of Oq, . . . , a^-x and will thus be 
known quantities for the graph, which one can correct for with appropriate phase gates. 

The model we will consider is a pseudo distance-regular network consisting of N sites 
labeled by {1, 2, A?"} and diameter D. In Ref. [1], we introduced the PST of a qubit in 
terms of the SU{2) generators. Let us now consider a state with d levels. First, we prepare 
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the generators for SU {d) systems and thereby introduce the Hamiltonians for a qudit system. 
The generators of SU (d) group may be conveniently constructed by the elementary matrices 
of d dimension, {epq\p, q e {1, d}}. The elementary matrices are given by 

i,^pq)ij ^ip^jqi ^ ^^jj ^ d, Cp . Cpp- (3-31) 

which are matrices with one matrix element equal to unity and all others equal to zero. These 
matrices satisfy the commutation relation 

^pq^ ^rs\ ^sp^rq ^qr^ps' 

There are d{d — 1) traceless matrices 

^pq ^pq ^qpi 

= -(epg - egp); l<p<q<d (3-32) 
z 

which are the off-diagonal generators of the SU (d) group. The d — 1 additional traceless 
matrices 

2 ^ 

Hm ^ -j= {Y, - msm+i}; m = 1,2, 1 (3-33) 

^J2m{m + 1) k=i 

are the diagonal generators so that we obtain a total of — 1 generators. SU (2) generators 
are, for instance, given as ax — — 621 + ei2, o-y = A21 = —i{ei2 — 621) and az — Hi — ei—e2- 
We now assume that at time t = 0, the qudit in the first (input) site of the network is 
prepared in the state jV'm)- We wish to transfer the state to the A^th (output) site of the 
network with unit efficiency after a well-defined period of time. As regards in the above 
argument, we choose the standard basis Cj = i — 0,1, ...,d — 1 for an individual qudit, 
and assume that initially all particles are in the state |0); i.e., the network is in the state 
|0) = |OaOO...OOb). Then, we consider the dynamics of the system to be governed by the 
quantum-mechanical Hamiltonian 

Hg^IJ: JmPm{\Y.^i ■ Xj + Ka. - \E\{^)]I,.), (3-34) 
m=0 i~i 
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where, is a — 1 dimensional vector with generators of SU {d) as its components acting 
on the one-site Hilbert space Hi, Jm is the couphng strength between the reference site 1 and 
all of the sites belonging to the m-th stratum with respect to 1, and PmS are polynomials 
given in (2-15) which are obtained using three term recursion relations (2-21) and the fact that 

— -j=====Q^. As it is seen from the Eq. (3-34), the terms of the hamiltonian for m > 1 
are nonlinear functions of \ " ^j- 

In the following we note that the term Hij :— Aj • \ in the hamiltonian (3-34), restricted 
to the one particle subspace (the subspace of the full Hilbert space spanned by the states with 
only one site excited), is related to the adjacency matrix of the corresponding graph. To do 
so, we write ifjj as follows 

iiii= E (Apf)®A+J^) + A;W^A;,(^))+|:i/«®i/^). (3-35) 

l<p<<j'<d m=l 

Before we proceed, one should notice that we have 

Cpg ® Crs = e(^p_i)a+r,(q~l)d+s (3-36) 

Then, from the fact that 

A+» ® A+(^) + A;» ® A,-(^) = 2(e« ® + eg ^ ej^)), 
and using the notation 

(m, n) = m -|- (n — l)d, 

one can obtain 

E (Ap+« ® A+(^) + A;» ® A;J^)) = 2 J2 h,,,),M + e^,,u.J (3-37) 
^<p<g<d '^<p<q<d 

Now, we evaluate the term Y^'^li Hm ® in (3-35) in terms of the elementary matrices 

Cpq as follows: First we note that 

2 fn rn m 

Hm®Hm^ —, — Ty{E E ^(p',p) " E [^(m+l.p) + e(p,m+l)] + e{rn+l,m+l)} ■ (3-38) 

mym + i) p=ip'=i p=i 
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The Eq. (3-38) can be rewritten as follows 

2 ™ 2m 2 

® = ^m(m + l) + ^^'^-^''-^'^^ + ^m(m + l)^_^,5>^^, 

Q m 

7 Z![e(m+i,p) + e(p,m+i)]}- (3-39) 



Then, one can show that 



m + 



d-l d 



Y,Hm®Hm = 2Y, e(p,p) - §/, (3-40) 

m=l p=l ^ 



for proof see Appendix A. 

Therefore, by using (3-37) and (3-40), H^j in (3-35) is written as follows 

d 2^ 

Hij = 2 ^ [e(g,p),(p,q) + e(p,5),(g,p)] + 2 ^ e(p,p) - -/. (3-41) 

l<P<9<<i P=l 

Now, one should notice that the permutation matrix which permutes two qudits, in terms 
of the elementary basis Cpq can be written as 

d d d 

P ^ Y Hp,q),{q,p) = H^{p,p)+ + (3-42) 

p,q=i p=l ^<P<q<d 

Then, the Eq.(3-41) takes the following form 

= 2Py (8) IaN-2 - ^ V , (3-43) 

where, denotes the permutation operator which permutes i-th and j-th sites and I^n is 
X identity matrix, where N is the number of vertices or sites {N :— \V\). 
Now, we denote a state in which the i-th site has been exited to the level u by = 

|0 O^j^O ... 0). Then, the permutation operators Pij only changes this state so that the 

i 

number and type of excited local states is fixed. This implies that the hamiltonian Hq can be 
diagonalized in each subspace S^''^ spanned by the vectors i = 1, . . . , A^, for z/ = 1, d — 1. 

We will refer to the states with only one site excited as one particle states and the subspace 
spanned by these vectors comprise the one-particle sector of the full Hilbert space. Then, the 
whole one particle subspace S can be written as 

5 = 5(l)©5(2)©...©5(<i-l). 
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In the other words, in d'^ dimensional Hilbert space, we deal with d — 1 one particle subspaces 
(recall that, each of these subspaces has dimension N). In the case of PST of a qubit, we have 
only one one-particle subspace of dimension N. 

One can easily show that, the operator J2i^j Pij restricted to the one particle subspace S^^\ 
can be related to the adjacency matrix A, as follows 

N 

E Pij - Ef'^max - i^imi"^ = ^ + (1^1 - «max)/iv, (3-44) 

where Kmax max{«;(i), i — 1,2, ... , N}, \E\ is the number of the edges of the graph, and 
H^^^ is the projection operator 7 (g) ... (8) 7 (8) Hi ®7 . . . 7 with 77; defined as in (3-33). For 

i 

regular graphs, where we have = k lor all i = 1,2, . . . , N , the equation (3-44) reads as 

TV — 2 

^P,^=A + k{^—)I^, (3-45) 



in which we have substituted \E\ = 

Then, by using (3-43) and (3-44), the hamiltonian in (3-34) restricted to each one particle 
subspace S^'^\ for u — 1,2, . . . ,d — 1, can be written in terms of the adjacency matrix A, as 
follows 

Hg = IJ2 JmPm{Y.P^3 ' + " |^|(^)]7^) = J^P^(A). (3-46) 

^ m=0 ^ ^ m=0 

For the purpose of the perfect transference of a qudit, we consider pseudo-distance- regular 
graphs with kd = \^d{o) \ = 1, i.e., the last stratum of the graph contains only one site. Then, 
we impose the constraints that the amplitudes {4>i^^\e~^^^°\(j)'Q^) be zero for alH = 0, 1, 7) — 1 
and v — 1,2, . . . ,d — 1, and (</i'^^|e^*^*°|0o'^^) = e*^, where 6 is an arbitrary phase. Recall that 
we have 

|#) = |^o...o), |#))=p,(A)|#); ^ = i,2,...,d-i. 

Therefore, the amplitudes {(f)\'^^\e~^^^°\(f)Q'^), lor i = 0,1, . . . , D must be evaluated. One can 
easily show that these amplitudes are independent of the value of v, i.e., it suffices to evaluate 
them for one choice of u, say v — 1. 
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As it has been shown in Ref. [1], the above mentioned constraints (by choosing 1^=1, the 
other choices give the same results), are equivalent to the following ones: 



fe=0 

Denoting e~^**" ^m=o by f]k, the above constraints are rewritten as follows 

D 

^ Pi(a;fc)%7fc = 0, i = 0, 1, - 1, 



A;=0 



D 



k=0 



(3-47) 



The Eq.(3-47) can be written as 



?7o7o 
\ VdId j 



( n \ 



(3-48) 



with Pij = Pi{xj). From the orthogonality relation (2-20) and using (2-25) , one can obtain 

Si, = {<p,\<Pj) = J2P^i^lHPA^l) ^ pwp' = i ^ pw'/'{pw'/y = i ^ 

(PW^V2)-1 = (piyl/2)t ^ p-1 = Wpt^ 

with W :— diag{'yo: 7i, ■ ■ ■ , 7d)- Therefore, P is invertible and the Eq. (3-48) can be rewritten 
as 



Volo 



p-1 



f ^ 





Ad 



^ ^ 



WP^ 



(3-49) 



The above equation implies that rjk'jk for k — 0,1, ...,D are the same as the entries in the 
last column of the matrix P~^ — WP* multiphed with the phase e*^, i.e., for the purpose of 
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PST, the following equations must be satisfied 

r7fe7fe = 7fce"''*°^-°^'"'''"^"'=^ = e*'(^P*)fcD , ior k = 0,1, D. (3-50) 
One should notice that, the Eq.(3-50) can be rewritten as 

(Jo, Ji,..., Jd) = --^[^+(2/o+/(0))7r, ^+(2/i+/(l))7r, . . . , e+{2lD+ f{D))T:]{WF'), (3-51) 
or 

Jk^-^Y.[0 + C^lm + f{m))7r]{WPXk, (3-52) 

where ior k — 0,1, . . . , D are integers and f{k) is equal to or 1 (we have used the fact 
that 7fc and (VFP*)^.^ are real ior k — 0,1, . . . , D, and so we have 7^ = K^P*)?;^!)- result 
(3-52) gives an explicit formula for suitable coupling constants so that PST between the first 
node (|0o)) and the opposite one (l^i?)) can be achieved. 

4 Examples of pseudo-distance-regular networks 

1. The networks Gn 

The networks Gn presented in [28], consist of two balanced binary trees of height n with the 
2" leaves of the left tree identified with the 2" leaves of the right tree in the simple way shown 
in Fig. 1 (for n — 2). The number of vertices in G„ is 2"+^ -|- 2" — 2. For the purpose of PST 
over Gn, we prepare the initial state to be transferred, at the left root of the graph and wants 
to calculate the suitable strength coupling constants so that the probability of the presence of 
the initial state at the right root be equal to one for some finite time tg. One can show that Gn 
is a pseudo-distance-regular graph with {2n + 1) strata, where stratum j consists of 2^ vertices 
for j — 1, 2, n + 1 and 2^^"'+^~'') for j — n + 1, 2n + 1. Therefore, its QD parameters are 

ai — 0, i — 0,1, . . . ,2n; uJi — 2, i — 1,2, . . . ,2n. 
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Then one can show that [19, 20] the polynomials Pi{x) are given by 

1 



Piix) 



=Qi{x) = Ui 



X 



(4-53) 



where UiS are the Tchebishef polynomial of the second kind. Then, xi^s (the roots of the 
Tchebishef polynomial Q2n+i(x) — •\/2^"+^C/2n+i(^)) and the coefficients 7/ are given by 



cos 



(/ + l)7r 



7/ 



2(n+ 1)' 

(-1)'+! . {1 + 1)tt . (2n+l)(/ + l)7r 



sm ■ 



sm 



, / = 0, l,...,2n. 



n + 1 2(n + l) 2(n + l) 

Therefore, we have Pij — Ui{-^) — [/^(cos or equivalently 



(4-54) 



Ui{^) ^i(cos^) ... t/i(cos(|g^) 



Furthermore, as it is seen from (4-55), the matrix P is a polynomial transformation [29]. By 
using the result (3-52), we obtain 

1 ^ (-1)™+! . (m + l)7r . (2n + l)(m + l)7r,, , (m + l)7r, 

2^0 ^0 ^+1 2(n-Fl) 2(n-M) 2(n-hl) 

(4-56) 



(4-55) 



for A; = 0, l,...,2n. 

In the following we consider the case n = 2 in detail. In this case, we have 



Xq 



\/6, xi = \/2, X2 = 0, X3 = -\/2, 



X4 



7o = 74 = 7i = 73 = 72 = g- 
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Then by using the equations (4-55) and (3-49), we have 



/ 1 


1 


1 


1 


1 \ 
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a/3 

V "J 


1 \ 
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-1 
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3 





-3 


-3 


2 





-1 





2 


12 


4 





-4 





4 


V3 


-1 





1 






3 


-3 





3 


-3 


. 1 


-1 


1 


-1 


1 J 




.1 


-V3 


2 




1 J 



Now, using the Eq. (4-56), one can obtain the following suitable coupling constants in order 
to PST between the most left and the most right nodes be achieved: 

_ e + 2n/3 J _ J _ 71 J _ J _ 71 

2. The modified G^-type networks 

This type of networks consist of two balanced binary trees of height n + 1. For the graphs G„, 
set entrance and exit to the left root and the right root of the trees, respectively (see Fig. 2 
for n = 2). These networks are pseudo distance regular networks which are characterized by 
the following intersection array 

{bo, 6i, . . . , 62n-i; ci, C2, . . . , C2n} = {3, 2, 2, . . . , 2, 1, 1, . . . , 1; 1, 1, . . . , 1, 2, 2, . . . , 2, 3}. 

^^^^^■V^^^^^ ^^^^^"V^^^^^ ^^^^^"V^^^^^ ^^^^^"V^^^^^ 

n—1 n n n—1 

Then by using (2-8) and (2-12), one can obtain 

K = Ki — 3, Ki — K2n—i = 3x2* ^, i = 1, 2, . . . , 71 K2n — Ij 

ctj = 0, = 1; a;i = a;2n = 3, a;^ = 2, i = 2, . . . , 2n - 1. 

Now, by using the QD parameters Ofj and cUi, and the recursion relations (2-21) and (2-23), 
one can obtain the Stieltjes function for any given n and spectral distribution ii{x) for the 
modified G„-type network. In the following we consider the case n = 2 in detail. 
In this case, we have 

ai — 0, i — 0, 1, 2, 4, as =1; uji = cu^ = 3, 002 = ^-^3 — 2. 
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Then by using the recursion relations (2-21) and (2-23), one can obtain the stieltjes function 
as 

x'^ — Ix^ — + 2x + Q 
" x(x-3)(x3 + 2x2-4x-7)' 

By using the recursion relations (2-21), we obtain 

Qi(x) = X, Q2{x) = x^-S, Qz{x) = x{x^-b), Qi{x) = x^-x^-TxVSx+G, ^5(0;) = x{x'^-x^-lQx^+bx+21). 

(4-57) 

Then, from the fact that, x;, / = 0, 1, 2, 3, 4 are roots of Q^{x) (the denominator of G^[x)), we 
obtain 

xo = 0, xi = 3, ~ -2.4728, x^ ~ -1.4626, x^ ~ 1.9354; 
70 = ^, 71 = ^, 72 - 0.3101, 73 - 0.1786, 74 ~ 0.1872. 
By using the Eq.(2-21) and the fact that Py = Pi{xj) = , / Q.^ixj). we obtain 

'ill 1 1 

^/3 -1.4277 -0.8445 1.1174 

-1.3417 1.3930 -0.3850 0.1492 

2V3 -0.4594 0.6974 0.4046 

1 2 0.5572 -1.4304 -0.6270 

^ 0.2869 -0.0155 -0.3725 -0.1982 0.2133 ^ 

-0.0303 -0.0157 0.0493 0.2495 0.0965 

0.2509 -0.1527 0.2961 -0.0807 0.1465 

0.0712 -0.3078 -0.1565 0.2957 -0.2813 

0.4213 0.4916 0.1836 -0.2664 -0.1750 
Now, using the Eq. (3-52), one can obtain the following suitable coupling constants in order 
to PST between the most left and the most right nodes be achieved: 



V 



Jo 



Q + 0.49 257r 



0.18387r 
2tn 



0.027l7r 
2tn 



0.02937r _ 0.45 637r 
-, J4 - 



2tn 



2to 
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3. Pseudo-distance-regular network derived from Icosahedron network 

This network has 12 nodes with 

{bo, 6i, 62; ci, C2, C3} = {5, 2, 1; 1, 2, 5}. 
Then by using (2-8) and (2-12), one can obtain 

K = 5, = 5, = 1, 

cco = 0, ai — a2 — 2, 0:3 = 0; a;i = 5, 002 = 4, a;3 = 5. 

The stratification with respect to the initial node 1 and the final node 12 produces the same 
strata (see Fig. 3). Therefore, the PST between the nodes 1 and 12 can be achieved. To do 
so, by using the recursion relations (2-21), we obtain 

Qi{x) = X, Q2{x) = x^-2x-b, Qsix) = x^-4x^-5x+10, Q^ix) = a;^-4a;^-10a;^+20a;+25. 



(4-58) 



Then, from the fact that, xi, / = 0, 1, 2, 3 are roots of Q4, we obtain 

^^51 1 
xo = -1, xi = 5, X2 = V5, ^3 = -V5; 70 = — , 71 = — , 72 = 73 = ^• 



By using the Eq.(2-21) and the fact that Py = Pi{xj) 



^ 1 111^ 

-75 ^ -1 1 
1 1-1-1 



(P)-i = WP* = 



12 



, / Q.Axj), we obtain 



^ 5 -^fl -\/5 5 ^ 

1 V5 1 

3 3 -3-3 
3-3 3-3 



Now, using the Eq. (3-52), we obtain the following couphng constants 
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4. 3- simplex fractal with decimation number 6 = 2 

Consider the (n — l)-simplex fractal with decimation number b — 2 [30]-[32] such that all of the 
n — 1 vertices (20. ..0), (020. ..0), (0...02) are connected to each other (see Fig. 4 for n = 4 



n-l 



n-l 



n-l 



). Then, the number of vertices of {n — l)-simplex fractal is N — C""^ + Ca""^ = n{n-i) ^^^j^ 
that the degree of each vertex is /t = 2(n — 2). Also, it can be easily shown that the network 
has 3 strata with respect to the reference node (200. ..0). For n — 4, the last stratum contains 
only one nodes so that the PST can be considered; In this case, the intersection array, size of 
strata and the QD parameters are given by 

{6o,6i;ci,C2} = {4,1; 1,1, 4}, 

K = 4, K2 = l, 

ao = 0, tti = 2, a2 = 0; lvi = 002 = 4. 
Then, by using the recursion relations (2-21), we obtain 



Qi{x)^x, Q2{x) ^ x"^ - 2x - 4, Qsix) ^ x^ -2x'^ -8x. 



(4-59) 



Then, from the fact that, xi, Z = 0, 1, 2, 3 are roots of Q3, we obtain 

111 

xo = 0, xi = 4, X2 = -2; 7o = -, 71 = -, 72 = -. 



By using the Eq.(2-21) and the fact that Pij = Pi{x 



P = 



V 



2-1 
-111 



(P)-i = WP* 



( 

1 



6 



V 



==Qi{x^\ we obtain 

3 0-3^ 
12 1 
2 -2 2 y 



By using the Eq. (3-52), we obtain the following constants in order to PST be achieved: 

2^ + 5a/3 



Jo 



1 J\ — 1 Ji — • 

4^0 3^0 12^0 
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5. Pseudo-distance-regular network derived from Hadamard network with 16 
nodes 

Consider the pseudo-distance-regular network given in Fig. 5. This network is obtained from 
the Hadamard network with intersection array {4, 3, 2, 1; 1, 2, 3, 4}. As Fig. 5 shows, the net- 
work is symmetric with respect to the initial and final (horizontal) nodes 1 and 16 e r4(l) 
and also with respect to the initial and final (vertical) nodes 6 and 11 G r4(6). One should 
notice that stratification of the network with respect to the nodes 1 and 16 produces the same 
strata. For stratification with respect to the node 1 or 16, we have 

K — A, K2 — 6, Ks — 4, K4 — 1, 

ai — 0, i — 0,...,A; cui — A, a;2 = 6, CU3 — 6, cu^ — A. 
(see Eq. (2-11)). Then, by using the recursion relations (2-21), one can obtain 

Qi{x) = X, Q2{x) = x'^-4, Qaix) = x^-lOx, Q^ix) = a;^-16a;^+24, Q^ix) = a;^-20a;^+64a;. 

(4-60) 

Then, from the fact that, xi, / = 0, 1, 2, 3, 4 are roots of Q^, we obtain 

3 1 1 

xo^O, xi^2, X2 = -2, 2:3 = 4, X4 = -4; 70 = 71 = 72 = 7, 73 = 74 = 7^- 

8 4 16 

Then by using the equations (4-55) and (3-49), we have 
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Now, using the Eq. (3-52), we obtain the following constants: 

_ e + bn/S J _ J _ TT 37r t Stt 
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PST of an arbitrary d-level quantum state over antipodes of pseudo distance-regular networks, 
was investigated. By using the spectral analysis techniques and algebraic structures of pseudo 
distance-regular graphs an explicit formula for coupling constants in the Hamiltonians was 
given so that state of a particular qudit initially encoded on one site can be evolved freely to 
the opposite site without any dynamical control. 
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Appendix A 

Proof of the Eq.(3-40): 

By using the Eq. (3-39), we have 

d—l d—1 2 2,711 '^^^ 



m=l m=l V / p=l m=l v / p,p'=r,p^p' 

2 

— — T ^hm+l,p) + e(p,m+l)]}- (A-i) 

We evaluate the first sum in the above equation, the second one can be evaluated similarly. 

^\ 2 ^ 2m ^ ^ d-i 2 ^ 2m 

i',^m(m+l)^^/(-'-) + - m(m + 1) + ^^^'^^^''-^'^ 

_^ 1 ' 2 4 2{d-l) _ 
- 1^ Cip-i)d+p + 3 + . . . + _ 2^ e(p_i)d+p + e<i+2 + ^ 62^+3 + ■ ■ ■ H ^ e<i — 

d 2fQ; 1) 1 1 ^ 

E ^ + 2(- - -)]e^a-i)d+a = 2(1 - 1/d) (ei + 6^+2 + . . . + 6^2) = 2(1 - 1/d) E e(p-i)d+p, 

a=i a a p^-y 

where, we have used the identity , = E*^"^ (- zr) = 2(i - i). 

' ^m—a m(m+l) ^m=a^m m+1' d' 

The second sum in (A-i) can be evaluated as 

d— 1 2 2 2 

^ ^^r^ , i\ E e(p',P) ~ ^ , 1 Ei^(™+l,P) + e(p,m+l)]} = E (^{p-l)d+p'- 

Therefore, we obtain 

d~l d t2 d d t2 

E ® = 2 E ^ip-i)d+p ~ E ^ip-i)d+p' = 2 E - -jl- 

m=l p=l ^ P,p'=l P=l 
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Figure Captions 

Figure- 1: Shows the network G2- 

Figure-2: shows the modified G2-type network. 

Figure-3: Shows the pseudo-distance-regular network derived from Icosahedron network. 
Figure-4: Shows the 3- simplex fractal with decimation number b — 2. 
Figure-5: Shows the pseudo-distance-regular network derived from Hadamard network 
with 16 nodes. 



002 



11 



15 



